Exact time autocorrelation function of the iV-spin classical Heisenberg equivalent 

neighbor model 
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We reduce the autocorrelation function Cn(t) of the equivalent neighbor model of N classical 
spins exhibiting Heisenberg dynamics and exchange coupling J to quadrature. As the temperature 
T — > oo, Cn(t) oc t~ N for Jt >> 1. At low T, the antiferromagnetic Cii(t) is a simple function of 
(JT) 1/,2 t, exhibiting strong frustration, but the ferromagnetic Cn(t) oscillates in a single mode the 
frequency of which approaches NJ as T — > 0. We conjecture that as T — > oo, the near-neighbor 
correlation functions of TV-spin classical Heisenberg rings are simply obtained from these results. 

PACS numbers: 05.20.-y, 75.10.Hk, 75.75.+a, 05.45.-a 



Recently, there has been a considerable interest in 
the physics of magnetic molecules. Q These consist of 
small clusters of magnetic ions imbedded within a non- 
magnetic ligand group, which may crystalize into large, 
well-ordered single crystals of sufficient quality for neu- 
tron scattering studies. Each type of magnetic molecule 
is characterized not only by the number N of spins in 
the molecule, but also by the spatial configuration of the 
spins. Many examples of rings exist, |2l but there 
are also examples of denser clusters. pOl Usually, the 
spins interact mainly via the Heisenberg exchange inter- 
action. Although interest in the dynamics of Heisenberg 
spin rings has been strong for many years, most of the 
work involved numerical simulation of the two-spin corre- 
lation functions. Recently, however, exact results for the 
Heisenberg dimer, isosceles and equilateral triangle, and 
four-spin ring were presented, (|, |[ but larger rings 
can not be solved using this technique. 

The equivalent neighbor model is the simplest model 
for the dynamics of a nanomagnet, first proposed by 
Kittel and Shore. || In this model, every spin interacts 
equally with all of the others. Although the Heisenberg 
dynamics of the N —* oo limit at infinite temperature 
T were obtained previously, they have not yet been 
solved for 5 < N < oo. Although this model is over- 
simplified for direct comparision with experiment when 
A > 5, its exact solution could serve as a benchmark for 
future theoretical and experimental studies. 

In this letter, we present an exact single integral rep- 
resentation of the autocorrelation function Cn(t) for the 
equivalent neighbor model with arbitrary A and T, both 
for ferromagnetic (FM) and antiferromagnetic (AFM) ex- 
change couplings. As T — > oo, an accurate asymptotic 
1/N expansion is found. As T -> 0, the AFM C n (t) 
reduces exactly to a simple function of (JT) 1 / 2 ^ for arbi- 
trary A. For the FM case at low T, Cn(t) oscillates in 
a single mode with a frequency that approaches NJ as 
T — ► 0, and a shape that fits a scaling relation. From ad- 
ditional studies of three- and four-spin rings with various 
amounts of bridging spins that interact only with the ring 



spins, we conjecture that the long-time T — > oo asymp- 
totic limit values of the near-neighbor correlation func- 
tions for equivalent neighbor coupling and near-neighbor 
coupling in classical Heisenberg A-spin rings are identi- 
cal. 

The Hamiltonian of the A-spin classical Heisenberg 
equivalent neighbor model is H = —Js 2 /2, where the 
total spin s = X^i s «W 1S a constant. The dynamics 
are given by dsi(t)/dt = Jsi x s, where Sj = |sj| = 1. 
Si(t) may be written as Si(t) = Si us + Si i j_[xcos(st*) — 
ysin(st*)], where t* = Jt, s = x x y, si ^(s,x) = (s 2 + 
I — x 2 )/(2s), x = |s — Si|, and s\ ^ + s\ ± = 1. Letting 
a = J / {2ksT), where fc# is Boltzmann's constant, the 
autocorrelation function C\\{t) = (si(t) ■ Si(0)) may be 
written 
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Cii(t) = — V N -i{x)dx I st 

Z N Jo J\x-1 



■ ts 'ds x 

x[s 2 ; ||(s,x) + cos(st*)s 2 ± (s,x)], (I) 

where the partition function Zn = J Q 2sT>N(s)e as ds 
and T>n(s) is the A-spin density of states. T>n{x) is 
piecewise continuous, 

E[(N-l)/2] 

V N {x) = B(x) ®(N-2p-x)x 

x6(i- N + 2p + 2)d N ^ 2p (x), (2) 
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where Q(x) is the Heaviside step function, (^) are bi- 
nomial distribution coefficients, and E(x) is the largest 
integer in x. The near-neighbor correlation function 
Ci2(t) = (si(t) • S2(0)) is simply obtained from Cn(i) 
and the sum rule (s 2 )/A = I = C n (t) + (A - l)C 12 (t). 

It is useful to write 5Cu(t) = Cn(t) — limi_ 00 Cn(t). 
Since Cn(0) = 1, it suffices to obtain 6Cu(t) for all t. Af- 



2 



ter some algebra, we have reduced SCn(t) to quadrature, 
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To obtain the long-time asymptotic behavior for \ot\ <C 
1, we integrate Eq. (Eh by parts N — 1 times. We find 
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(11) 



For N — M 2 , one of the terms in Eq. ( |l0| ) vanishes, as 
for the four-spin ring, [p] More important, the leading 
long-time behavior oc 1/t* , which falls off very rapidly 
for large AT. Although Eq. (ffo|) fails for \a\ > 1, its 
behavior for a ^ is useful to illuminate the dramatic 
difference between the long-time asymptotic behavior for 
the FM and AFM cases. For the FM case, a > 0, the 
dominant behavior oc cos[W(i* — 7r/2)], as all of the spins 
oscillate together. In the AFM case, a < 0, the dominant 
long-time asymptotic behavior is given by the smallest 
possible oscillation frequency. As shown in the following, 
Cn(t) is always finite as T — > 0. 

As T, N — > oo, Liu and Miiller obtained an analytic 
form for Cn, but their procedure was not extendable to 
finite N. |9) From the integral representation, T>n(x) = 
J °° 2p x ~ N sin N p sin(px) dp/ n, we expand sin N p/p N as 
exp[-A^(p 2 /6 +p 4 /180 +p 6 /2835 + ...)] to obtain a 1 /N 
asymptotic expansion for C\\{t) as T — ► oo, 

lim Cn(t) - [l + 2e- p (l-2i 2 )]/3 
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where t = Jt[(N - l)/6] 1 / 2 . From Eq. @, Cn(0) = 1 
through order 1/N 2 , as required, and Cn(t) has the cor- 
rect A" — > oo limit. [0] However, Eq. (|lj) is highly accu- 
rate for intermediate Af values as well. In Fig. 1, we have 
plotted limT^oo C n (*) for AT = 5, 6, 7, 8, 9, 25, 50, 100, oo 
as a function of Jt[(iV - l)/^] 1 ' 2 . The curves for AT = 
25, 50, 100, oo were obtained from Eq. (p^), whereas the 
other curves were generated from the exact solution, Eqs. 
(Bj)-(^). Equation (|l2|) is very accurate at short and 
long times, even for N = 5. We note that for t < 1, 
liniT^oo Cn(t) is a universal function of Jt^jN — 1, as 
suggested by the 1 /N expansion. Although the N — > oo 
curve is exact, and the curves for N = 50, 100 are quite 
accurate, the curve for A^ = 25 exhibits a small (« 1%) 
inaccuracy for 1.2 < t < 2.2. 

At arbitrary AT, i and AT 3> 1, f — » oo, respectively, 



lim Cn(t) = 

T — >oo 



l/N+(N-l)[6 N + f N (t)},{13) 
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In the inset of Fig. 1, we have compared the exact values 
of Sn with those obtained from Eq. ([To] ) for 3 < A" < 65. 
For A" = 3, the differences are noticeable on this curve. 
As N increases from 5 to 65, the relative difference de- 
creases rapidly from 0.14% to 0.2 ppm. Since the exact 
formula for 8jy with N odd (even) contains non-trivial 
contributions of logarithms of all prime numbers from 3 
to A" (2 to N/2), respectively, the extremely good nu- 
merical agreement between the exact 5n and Eq. (^) is 
remarkable. 
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FIG. 1: Plots of limT^ooCii(t) versus \J\t[(N - 1)/6] 1/2 , 
for TV = 5,6,7,8,9 from Eqs. (4)-(9) (solid), and for N = 
25, 50, 100, oo from Eq. (12) (dashed). Inset: Plot of the 
exact (x) and asymptotic values, Eq. (15), of Sn- 
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FIG. 2: Plots of limt^oo Cu (t) versus \a\ for N = 5,7, 10 and 
for both ferromagnetic and antiferromagnetic couplings. 



In Fig. 2, we plot the long-time asymptotic lim- 
its lim t ^oo C\\(t) for both the FM and AFM cases, for 
N = 5,7,10, as functions of \a\. The FM values ap- 
proach 1 as T — > 0. But, the low-T AFM values rapidly 
approach their T = asymptotic limit, |. This indicates 
the spins are strongly frustrated, since one might naively 
expect the AFM lirnt_, 00 Cn(t) — 0, as for the four-spin 
ring. H In that case, the unfrustrated alternating spin 
direction configuration was possible. In this AFM case, 
however, the N > 3 spins are as frustrated as T — > as 
are an infinite number of them as T — ► oo. 

We now examine the low-T AFM dynamics. In this 
case, the dominant contibutions to the numerator of Cu 
and to Zjq arise from s <C 1. The leading contributions 
to these integrals are P/vs 2 e _ ' Q ' s and QArs 2 e~' Q ' s , re- 
spectively, where P/v and Qn are complicated functions 
of N, but for arbitrary N > 3, Pn/Qn = §• Thus, as 
T — ► 0, the Fourier transform 6Ch(li) of 5Cn(t) attains 
the exact uniform asymptotic form, 



3^ 



o) 2 e 



(16) 



where w = |a| -1 / 2 a;/| J|. For A > 5, corrections to Eq. 
( |l6| ) are of C(|a| _1 ), whereas for N = 4, they are of 
0(|a;| _1 / 2 ). In Fig. 3, we show the frequency dependence 
of the low-T AFM mode. Since C n (t) for fixed A is a 
uniform function of |J|i/|cn| 1//2 , the amplitude of SCh(uj) 
must also be scaled by |a| -1 / 2 for the curves to nearly 
coincide. Here we plot the curves at a = —5,-10 for 
N = 4, 7, 10, obtained from our exact solution, Eqs. (Q) - 
(Q). The exact asymptotic form, Eq. ( |l6| ) is also plotted 
for comparison. The N = 4 curves nearly coincide with 
the others, but do show some small dependence upon T. 
However, the N = 7 and 10 curves are nearly identical 
and independent of T, and are nearly indistinguishable 
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FIG. 3: Plots of the low-T \a\- 1/2 SCu(iu 
for the antiferromagnetic cases with JV = 4, 7, 10 at a 
—5, —10. Also shown is the exact T — *• limit, Eq. (16). 



u/\J\ 



from the asymptotic curve for \a\ — 10, when scaled in 
this manner. 

Inverting the Fourier transform, the low-T AFM Cu(t) 
asymptotically approaches 



Cn(t) T :i [l + 2e- 4 (l-2F)]/3, 



(17) 



where t = t* /(2\a\ 1 / 2 ). Hence the AFM lim T ^oCn(t) 
is a uniform function of (JT) 1 / 2 f, independent of N. 
Moreover, it has precisely the same form as does 
limjv^—,00 Cn(i), Eq. (p"2|), pictured in Fig. 1, except 
for the different scaling factor. 

We now turn to the FM case as T — > 0. From Eq. 
(|l0|), we expect that all of the spins will be oscillating 
together with frequency NJ. But, this is strictly true 
only in the T — > limit. At finite a, the oscillation 
frequency deviates from this value slightly. To determine 
more precisely the actual nature of the mode, we note 
that for a > 1, the peaks in the integrand and in Zn both 
occur within A — 2 < s < A. AsT^O, we evaluate Zm 
by integration by parts, leading to the exact low-T limit 
of the FM Fourier transform, 



SCu{uj) 



A N f N (L0*)/f N (LU* N ), 



(18) 



where A N = 2e[(A - l)/e] w /A!, uj* = uj/J, and u* N , the 
position of the exact maximum in /jv(w*), is an eigen- 
value of 2as — (A — 2)/(A - s) + d\n[g (s)]/ds = 0. f N 
and go are given by Eqs. (||) and @, respectively. Ap- 
proximately, oj* n ps A — (A — l)/(2Aa), and fi^(u>*) has 
a width characterized by a normal Gaussian distribution 

Thus, 



parameter cr/v 



An exp[— (ui* 



co N r/(2a%)}. (19) 
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FIG. 4: Plots of the low-T SCn(cu)/A N versus (u/J - 
io* N )Na/VN - 1, where cj n = N — (N — l)/(2Na) and 
A N = 2e[(N - l)/e] N /N\, for the FM cases with N = 8, 
a = 5, 10, and N = 5, a = 5, 10, 28. 



In this Gaussian approximation, u)^, <jn, and A jy are cor- 
rectly given to their respective leading orders in 1/a, but 
the skewness of the peak is not accurately represented. 
In Fig. 4, we have therefore plotted the scaled low-T 
5C n (w)/A N versus (<j* - to%)Na/^/W^T, for N = 8 
at a = 5, 10, and for N = 5 at a — 5, 10, 28, using our 
exact formulas, Eqs. (Q) - (||). The A^ = 5 curves nearly 
coincide, as do the N — 8 curves, and the scaled height 
approaches 1 as T — > 0. As mentioned previously, the 
asymmetry of the curves decreases with increasing N, a 
feature not incorporated in the Gaussian approximation. 
Neverthess this scaling procedure works remarkably well, 
and N — 5 is close to the large A" limit, even for the FM 
case. 

Inverting the low-T FM Fourier transform, Eq. (|l9|), 



Cu(t) 



1 



Bn 
a 



1 _ e -t' 2 a 2 N /2 C0S( V^*)1 ; (20) 



where B N = A Ny /n(N - l)/2/N. As N -> oo, B N 
1/N. The deviation of the long-time limit of Eq. (fee 
agrees to within 1% with our results presented in Fig. 
2. The decay in time of the mode in this approximation 
is Gaussian, with a lifetime tjv = l/fx/v- Thus, as T — > 
0, the oscillations are characterized by a frequency that 
approaches NJ linearly in T, a lifetime that diverges as 
1/T, and an amplitude that vanishes as T. 

Finally, we remark that our exact results for 
limr—,00 Cn(t) in the AT-spin equivalent neighbor model 
may also lead to the correct T — > oo limit of the near- 
neighbor correlation function C\i (t) of Af-spin rings with 
classical Heisenberg near-neighbor exchange coupling. 



From Eq. (|13|) and the sum rule, we conjecture that, 
as for the A-spin equivalent neighbor model, 



lim C 12 (t) = 1/N - 6 N - f N (t) 



(21) 



for the A^-spin ring. In support of this conjecture, we 
have studied the A^-spin systems consisting of M spins 
s'j coupled to N — M spins Sj on a ring. Letting s = 

SiIo~ M Sj and s' = YljLi s 'j> tne Hamiltonian is 



N-M 



H 



J E 

i=l 



St+i - J's ■ s', 



(22) 



where sjy_Af + i = si, which is integrable. Q Setting 
C'(t) = (s'(i)-s'(O)), we found for N—M = 3, M = 1,2,3, 
and for N - M = 4, M = 1, 2, that 

lim C'{t) = M 2 /N + M(N - M)[5 N + f N {t)}.{23) 

T — >oo 

Most significant of these is the result for the four-spin 
ring with M = 2. We also showed previously that Eq. 
( pl| ) is valid for four- and three-spin rings. [|[ ^ 

In summary, we have solved exactly for the time auto- 
correlation function C\x(t) of a spin in the A^-spin equiv- 
alent neighbor model. As T — > oo, Cn(t) oc t~ N for 
Jt 3> 1, curves for different N exhibit very similar shapes 
when plotted as functions of | J\t[(N — l)] 1 ^ 2 , and an ac- 
curate 1/N asymptotic expansion is found. At low T, the 
antiferromagnetic Cn(t) is a universal function of tT 1 / 2 , 
independent of N, and exhibiting strong frustration. For 
the ferromagnetic case at low T, there is a single mode 
with a peak position at w NJ — ksT(N — 1)/N, and a 
shape that fits an accurate scaling relation. We conjec- 
ture that the T — ► oo limit of the near-neighbor correla- 
tion functions C\2(t) for the A^-spin equivalent neighbor 
model and for the A^-spin ring may be identical. 
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